'Spin-Spin' Hall Effect in Two-Dimensional Electron Systems with Spin-Orbit 

Interaction 



o 
o 

(N 
> 

o 
in 



Ankur Sensharma and Sudhansu S. Mandal 

Theoretical Physics Department, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India 

(Dated: February 2, 2008) 

We calculate spin-Hall conductivities in the two dimensional electron systems with Rashba spin- 
orbit interaction. The salient feature is that, apart from the usual spin-Hall conductivity a% y 
which corresponds to the induction of out-of-plane spin-current due to the application of transverse 
charge current, there is a novel spin-Hall conductivity a^r y which arises due to the induction of 
transverse spin-polarized current in the transverse direction by the application of in-plane spin- 
polarized current. This phenomenon which we call as 'spin-spin' Hall effect is a spin analog of 
conventional Hall effect, but with no magnetic field. This contribution may be understood through 
the spin-diffusive equation. 



PACS numbers: 72.25.Dc, 72.80.Ey, 72.25.Rb, 71.70.Ej 
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INTRODUCTION 

Following pioneering proposal of a device called spin 
field effect transistor (SFET) by Datta and Das[l], there 
has been growing interest in the field of spintronics[2] 
which is the science of coherent manipulation of spin, 
because of its potential applications in spin-memory and 
quantum computing devices [3]. This device consists of 
a two-dimensional electron system (2DES), controlled by 
gate voltage in a semiconductor heterostructure with spin 
polarized contacts. A spin entered from a spin-polarized 
source in the 2DES, precesses due to the Rashba spin- 
orbit (SO) interaction arising from the lack of structural 
inversion symmetry in the semiconductor heterostruc- 
tures. The spin of the electrons can be made to trans- 
port coherently and the spin-polarized drain whose spin 
polarization is parallel to the same as source can de- 
tect the transportation of spin. The spin current can 
be manipulated by controlling gate voltage since the 
Rashba SO coupling depends on it. However, the mo- 
mentum relaxation due to elastic scattering of the elec- 
trons from impurity potential leads to spin relaxation and 
thereby destroys spin coherence. Although SFET has 
not been achieved yet, the modulation of SO coupling by 
gate voltage has been observed in InGaAs/InAlAs and 
GaAs/AlGaAs heterostructures[4-6]. This motivates to 
study the system of electrons with SO interactions. 

A remarkable consequence of the SO interaction is the 
spin-Hall effect where an electric field induces a trans- 
verse out-of-plane spin-current and thereby a spin imbal- 
ance takes place between two opposite edges of the sam- 
ple. Ever since the prediction of intrinsic dissipationless 
spin-Hall effect in hole-doped semiconductor systems [7], 
this effect has drawn intense theoretical as well as exper- 
imental activities. Although Sinova et al [8] have pre- 
dicted universal (independent of the strength of Rashba 
spin-orbit interaction) dc spin-Hall conductivity (SHC) 
a xy — e/8-7r in a pure 2DES, the issue of SHC in the 
presence of non-magnetic impurities remains highly con- 



troversial. A number of analytical works [9-11] based 
on Kubo formula and quantum kinetic equations suggest 
that a xy vanishes for any amount of elastic disorder in 
the diffusive transport regime. This dramatic vanish- 
ing of SHC in the presence of disorder is argued to be 
a more generic phenomenon [12] as a z xy also vanishes for 
any amount of magnetic field in a pure system. On the 
contrary, experiments [15-17] seem to suggest that there 
is nonequilibrium spin accumulation in the edges, trans- 
verse to the direction of application of charge current. 
Numerical studies [13] based on the Landauer-Buttiker 
approach in mesoscopic systems with leads cannot shed 
any light about SHC in the thermodynamic limit because 
of the possibility of edge effects near the contacts. How- 
ever, a numerical calculation based on the Kubo formula 
in the lattice model [14] suggests that SHC is finite in 
the presence of disorder, while it vanishes in the ther- 
modynamic limit. Nomura et al [18] have numerically 
shown that the vanishing of intrinsic o~* y is peculiar to 
the linear-momentum Rashba model which has been con- 
sidered in the previous studies. They have shown that 
for Rashba model with cubic momentum which is the 
case for two dimensional hole gas, the intrinsic a xy is 
nonzero and consequently intrinsic edge-spin accumula- 
tion takes place in the systems of two dimensional hole 
gases [16, 19]. On the other hand, the edge-spin accu- 
mulation observed [15, 17] in the system of two dimen- 
sional electron gas where Rashba model with linear mo- 
mentum is important, is argued to be extrinsic by the 
'skew-scattering' mechanism [20-22]. While the issue of 
intrinsic versus extrinsic spin-Hall effect is not clear yet, 
we show below that the linear Rashba model provides a 
different but robust kind of spin-Hall effect in the two 
dimensional electron systems. 

In this paper, we calculate the spin-Hall conductivities 
in the charge-spin space using Kubo formula in section 
III. A charge current induces an out-of-plane component 
of spin current in the transverse direction. This usual 
ac spin-Hall conductivity a xy is found to be the same 
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as in a kinetic equation approach [10]. It however van- 
ishes in the dc limit [9-11]. We call this phenomenon as 
'charge-spin' Hall effect. Apart from this, we find that 
an in-plane spin-polarized current can induce a trans- 
verse spin-polarized current along transverse spatial di- 
rection. This phenomenon which we call 'spin-spin' Hall 
effect is a spin-analog of the conventional charge Hall ef- 
fect, but in the absence of magnetic field. The relevant 
Hall conductivity of this phenomenon a xy does not van- 
ish in the dc limit. We attribute this phenomenon to the 
diffusion of spin s y along y-direction due to a source caus- 
ing conventional s x spin-diffusion along x-direction. The 
corresponding diffusive equation is derived in section IV. 
The detailed derivation of diffusive propagator is avail- 
able in the Appendix A. In Section V, the spin-spin Hall 
conductivity has been argued to be robust. Relevance 
of other type of spin-spin Hall conductivities have also 
been discussed. The consequence of Dresselhaus SO cou- 
pling along with the Rashba coupling on spin-spin Hall 
conductivities has been discussed in Appendix B. 



HAMILTONIAN AND GREEN'S FUNCTIONS 

A system of two dimensional non-interacting electrons 
with Rashba spin-orbit interaction in the presence of non- 
magnetic impurities from which the electrons scatter elas- 
tically can be described by the Hamiltonian 



H = H Q + V(r): H 



P 

2m 



At) • p , 



(1) 



where p = — iV is the momentum operator, m is the 
effective band mass of electrons, A is the parameter for 
the strength of Rashba spin-orbit interaction, r) = n x a 
is a spin operator with <r being the Pauli matrices, n is 
the unit vector perpendicular to the plane of the system, 
and V(r) is the disorder potential. The eigenvalues of 
Hq are e k = k 2 /2m + sXk where k is the magnitude of 
the wavevector k of the electrons, s — ± is the index for 
spin-split subband s. The corresponding cigenstates are 
given by 



l 



( 



iXk/2 



^2 V se-^/ 2 



(2) 



where %k = n/2 — with 9 being the angle of k with 
x-axis. The Fermi momenta and the density of states 
at Fermi energy, ep, in the two subbands are k s F — 



\Jkp + m 2 X 2 — smX and V s = v(\ — s- 



:) respec- 



tively, where v = m/2-K is the density of states for each 
spin direction and fcjr = \J2mep is the Fermi momentum 
in the absence of SO interaction. The charge and spin 
currents are j Q = e a \ {a a , vk} where Vk = ^oo + At) is 
the velocity and &q is the unit matrix . Here a = 0-3 
represents charge and three spin directions respectively 



and eo = e (charge of the electrons), e\ = &2 = &z = 1/2 
(spin of the electrons). 

We next assume random disorder potential with the 
configurational average such that (V(r)) = and 
(V(r)V(r')) = 7<5(r - r'). The elastic life time of the 
electrons will then be r = l/(2n-fv). Therefore, the re- 
tarded (advanced) Green's function for the electrons with 
energy e and wavevector k can be expressed as 

9 (k ' £) ^L £ _ a±l /2r^ gs (k ' £) ' 



s=± 



s=± 



(3) 

where r)k = t) • k/k is the projection of the spin operator 
into the direction of k and £ k = e k — ep. 



HALL CONDUCTIVITIES 

The spin-Hall conductivities a* y , a^ y , and a\ y cor- 
respond to transverse induction of out-of-plane (z-axis) 
spin-current due to a charge current, transverse (in plane) 
spin-current due to a spin-polarized current, and paral- 
lel spin-current due to a spin-polarized current respec- 
tively. These conductivities can be obtained using Kubo 
formula: 



(4) 
(5) 
(6) 



where retarded off-diagonal current-current correlation 
functions can be written as 
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<M = / 7^Tr(i-(k)^(k,0) [%W + j%{u>)]g R (Kuj) 



dk 

(2V) 

( 7 ) 

Here the current j/3(k) corresponds to the bare vertex. 
This vertex gets corrected through Jp which is the self- 
consistent solution of the equation 

J/jM=7 / ^a' 4 (k',0)[j /3 (k') + J^M] e*(k». 

(8) 

It describes summation over infinite number of the ladder 
diagrams in a diagrammatic approach. 

To solve Eq. (8) and evaluate Eq. (7), we need to con- 
sider integrations, e.g., 



I 



dk 



#(k,0)#?(k,w) 



(2n) 2 ^ 

which we calculate using the approximation 

r 2ir 



f dk dO r s f°° , t 



d£k{l - 5 SS >) 

> 

(9) 
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i.e., the integrations are performed around momentum 
kp for the intraband contribution while the interband 
contributions are obtained through the expansion around 
kp and the variable of integration for the latter is £k = 
k 2 /2m — ep. This approximation is, however, neither bi- 
ased towards any of the subbands nor negligent of the 
existence of those. This attention is important towards 
our result. The self-consistent solutions are thus straight- 
forward and we find as 



JoM = 

JiH = 

J 2 (w) = 

J 3 (w) = 



-A 



26 2 



25 2 {1 - 2iur) - iujt{1 - iLUT) 2 
kp\ S 



2m ) 4<5 2 — iur{\ — itor) 
kp\ S 



Tj (10) 
(11) 



2m J 4<5 2 — iur{\ — itor) 
kp \ <5(1 — iuor) 



2m ) 2<5 2 (1 — 2iujt) — iwr(l — iujt) 



tr 3 £c 

(12) 
«tl3) 



for an arbitrary value of the parameter 8 = XkpT. While 
J and J 3 are due to both the intraband and the in- 
terband contributions, Ji and J 2 are due to interband 
contribution only. This means Ji and J 2 exist only be- 
cause of the formation of two subbands. In fact all of 
these vanish for A = 0, in consistence with the fact that 
the charge transport time and the momentum relaxation 
time are the same (no vertex correction) in a two dimen- 
sional system with short range scatterers. 

Using equations (4)-(13) we find the Hall conductivi- 
ties: 



IbJT 



2tt J 2S 2 (l - 2iLur) - iwr{\ - iujr) 2 
A \ A5 2 



(14) 



16ir J [A8 2 — iu>T(l — iu)t)][1 — itur] 
4 v {w) = 0. (16) 

The spin-Hall conductivity cr xy (u)) in Eq. (14) is identical 
to those derived from the quantum kinetic equation by 
Mishchcnko et a/[10] and in an approach of Kubo formula 
by Chalaev and Loss [11]. At ui = 0, the vertex correction 
term dramatically cancels [9-11] the bare contribution in 
a* y . The vanishing of a z y even in the weak localization 
regime [11] warns the presence of an exact property of the 
system under consideration. From a general argument, 
Chalaev and Loss [11] have shown a xy vanishes exactly 
for any finite amount of disorder. Since the spin is not 
conserved, the spin-Hall conductivity may not be directly 
related [12] to the observed [15-17] spin accumulation at 
the transverse edges. 

The most important result in this paper is the find- 
ing of the conductivity <J xy (uj). This suggests that an 
in-plane spin polarized current induces a transverse spin 
current with in-plane spin polarization but perpendicular 
to the former. If the source of the spin current is polar- 
ized along x-axis, then a detector, placed transverse to 



the direction of the applied current, with polarization 
along ±j/-axis can detect this novel current. Magnetic 
semiconductors may be used for spin injection as well as 
detection of the spin-polarized currents as described by 
Fiederling et al [23] in a light-emitting diode. Further 
techniques [15-17] developed for the detection of spin- 
accumulation in spin-Hall effect may also be useful for 
observing spin-spin-Hall effect. Unlike the conventional 
spin Hall conductivity a z xy , this new spin-spin Hall con- 
ductivity is finite at u> = and is given by 



(17) 



independent of the Rashba coupling parameter A. How- 
ever from the equation (15) we see that cr xy {uj = 0) 7^ 
only when A ^ 0. In other words, A — > before uu — > 
is logically the correct limiting procedure for dc a xy 
in the absence of spin-orbit interaction. The bare and 
vertex correction contribution to a^r„ in Eq. (17) are 



A 
16tt 



xy 



iS 2 



and _ W (t+I3^) respectively. Clearly 
the dominant contribution is due to the vertex correc- 
tion. The vanishing a\ y in Eq. (16) implies zero parallel 
spin-polarized current in the transverse direction. 



SPIN DIFFUSION 

We now derive the spin-diffusive equation for under- 
standing the above novel spin-spin Hall effect. The re- 
tarded charge and spin density correlation functions at 
frequency u> and momentum q can be written as 



Xa/s(q,w) = — 



dk 



2tt J (2tt) 2 



Tr a a g A (la,O)M^{k + q,uj)a 

(18) 

where a, (3 = 0-3 corresponding to charge and three 
spin directions, Gp(k, e) = apQ R (]s.,e)d-p, and Mp = 

Sjlo w i tn ™$ = ^0 an d 



m 



^^(k',0)mf 



(19) 

It is very easy to check that the Eq. (18) includes infi- 
nite number of ladder diagrams for the vertex correction. 
Although the expression of Xaf3 (18) appears to be un- 
common, we find it to be very convenient way to express 
for all orders of diagrams since we need to maintain the 
order of matrices inside the trace. The expression of Xa/3 
in (18) is formal. We however evaluate it for small q 
and lo such that ojt < 1 and q <C kp, and using the 
approximation in the expression (9) for the integration 
over electron momenta. We iteratively calculate the co- 
efficients of it's in ma's (19) and take sum of geometrical 
series of the coefficients. The readers interested in the 
details of the calculation may look into the Appendix A. 
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Expressing mfi* = J2l=o m i3a^ a and 
Xafi = 2ivujT{V afi - 5 a p) 



(20) 



in a standard form, we find the inverse diffusion propa- 
gator to be 



1 - m 00 -mu -m 2 2 

-m i l-m 10 -im 32 

~m 2 l-roao im Z i 

im\2 —im 2 \ 1 - m 30 



(21) 



where moo = 1 + r(iu> — Dq 2 ), mio = ni2o — moo ~ 
t/t|, to 30 = moo — T I T ^ w ith diffusion constant D = 
A/m, and = r(l + 4<5 2 )/4c5 2 and r| = r(l + 45 2 )/2<5 2 
being the out-of-plane and in-plane spin relaxation times 
respectively, tooi = mil = 2iq y S 3 /kp(l + AS 2 ), m 2 = 
m 2 2 = -(qx/q y )m i, m i2 = m 32 = -4:q x SA/k F (l+4<5 2 ) 2 
and m 2 i = TO31 = — (q y /q x )mi2- We thus find = 
-2q x q y AT/m(l + AS 2 ) 3 for w = 0. Clearly X*^ 1 and P^ 1 
are nonzero for small (5, contrary to the recent result of 
Burkov et al[24]. (See Appendix A for details). This 
is because the spin-relaxation provides an energy cut- 
off. These are indeed the responsible terms for driving 
transverse spin in the transverse direction as we shall 
show below. However, the other components of V^ 1 are 
in agreement with Rcf. 24. 

Transformation of the inverse diffusive propagator (21) 
into the real space and time leads to the transport equa- 
tions for spins with planer projection (neglecting gradi- 
ents of charge and out-of-plane component of spin densi- 
ties): 

f =(^-l) Sl - ?r ^ ¥W » + /,,22) 

f = ( DV2 - 4) S " - (TTW W " + (23) 

where S x (S y ) is spin density along x(y) direction and 
cc(y)-component of the spin current I x {Iy) is injected 
into the system. Although the coefficients of d x d v S y in 
Eq. (22) and d y d x S x in Eq. (23) appear to be nonzero for 
S = 0, these coefficients in fact vanish strictly at S = 0. 
The spin diffusive equations (22) and (23) are valid for 
wt <C AS 2 . Clearly, the injected spin current I x {I y ) not 
only transports the S x (S y ) component of spin along the 
direction of the current but also drives the S y (S x ) com- 
ponent of spin along the transverse direction. The finite 
dc spin-spin Hall conductivity (17) including sign arises 
as a manifestation of this fact. 



DISCUSSION AND SUMMARY 



however, have considered generally accepted definition of 
the spin-current operator, i.e., the symmetrization of the 
product of spin operator and the group velocity of elec- 
trons. Apart from this spin-current which depends on 
the translational motion of spin, there is a current asso- 
ciated with the rotational motion[25]. Albeit the defini- 
tion of current is not unique, the arbitrariness lies only 
in the terms proportional to the SO coupling parame- 
ter A. The kinetic parts[24] (which are proportional to 
the momentum of electrons and independent of A) of the 
spin-currents are, however, unique. The substantial con- 
tribution in a xy which arises due to vertex correction (17) 
is entirely due to the kinetic part of the spin currents, and 
it is not proportional to any power of A. We therefore 
believe that the 'spin-spin' Hall effect is robust in the 
system of two dimensional electron gas with the Rashba 
spin-orbit interaction. 

We have calculated a xy and a\ y which suggest spin- 
spin Hall conductivities when the induced spin-current 
will have respective spin-polarization perpendicular and 
parallel to the spin-polarization of applied spin-current. 
These spin-polarizations are in the plane of the 2DES. 
While the mechanism for usual spin-Hall effect in this sys- 
tem is mainly the momentum dependent spin-precession 
along with the spin relaxation, the mechanism behind 
spin-spin Hall effect is spin-diffusion. If the spin- 
polarization of the applied spin-current is out of the 
plane of the system, we find the corresponding compo- 
nents of off-diagonal current-current correlation functions 
(7) vanish: Ilg + 11^ - Uf y + Uf y = IE£ + Uf y - 
n 2 . 2 + U xy — 0. In other words, the spin-spin Hall effect 
is possible only when the applied spin-current's polariza- 
tion is in the plane of the system. 

In summary, we predict that an in-plane spin-polarized 
current can induce a transversed spin-polarized current 
along transverse direction. This phenomenon which we 
call 'spin-spin' Hall effect is a spin analog of conventional 
Hall effect, but with no magnetic field. The reason be- 
hind this phenomenon is the transverse spin diffusion 
along transverse direction due to the application of a 
source causing conventional spin diffusion. 
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DIFFUSIVE PROPAGATOR 



In this appendix, we present detailed calculation of the charge-spin density correlation function, and hence inverse 
diffusive propagator and compare with previous result [24] . The retarded charge-spin density correlation function can 
be represented by a sum of infinite number of ladder diagrams: 



1=0 

where (j) represents the order of the diagrams. Explicitly these are given by 



X$(«L") = ^Tr 



2tt 
iu) r 
2^ 



/ 7^2^ A ( k '°)^ fl ( k + < l' w ) 



dki 



(2tt) 2 y (2tt) 



a Q ^(k,0) 



x g A (k 1 ,0)a^ K (k 1 +q, W )^ K (k + q,^) 



(24) 



(25) 



(26) 



and so on. In the expression of Xafi, the contribution at u> — and q = for diagonal terms has been ignored. It 
is crucial to maintain the relative positions of the matrices inside the trace while performing the infinite series. To 
express all the orders in a compact form and to calculate those, we find it convenient to shift <jp to the extreme right; 
but in doing so, <7 K (k + q, uu) will be modified to (k + q, u) = &pQ R ()i + q, w)<7g. We thus find 



(0)/ \ 
Xa/j(q» w ) 



ILU 

2^ 
2^ 



Tr 
Tr 



_dk 

W) 

dk 



2 ^(k,0)^(k + q,w)^ 



r^(k,0) 7 



dkx 

(2^ 



g A (kx , 0)^(kx + q, w) } ^(k + q, w)^ 



(2)/ \ 



2^" 

7 



, o /^ (M) ( 7 /^ (k2 , 0) 



and so on. Defining 



m^(q, W ) = 7 y ^^(k'.O)™" A '(q, w)^(k + q, w) ; j > 1 



dk' 



Ci-i)/ 



(27) 



(28) 



(29) 



(30) 



with mjj ^ = <t , and summing over all orders (the lowest three orders given in Eqs. (27)-(29)), we rewrite the retarded 
charge-spin density correlation function as 



X«/3(q,w) = — Tr 



(31) 



Expressing 



= m/JO^O + m^iCTi + TO/32IT2 + m/33173 



(32) 



with /? varying from to 3, we obtain 16 components of mp a which are m o = 1 + r(ioj — Dq 2 ), m w = m 2 o 



TO °° — 1+4,52 ) m 30 — ™oo — 1^3525 m oi — 



2i 



</•< 



4<5A 



fc7(T+4P)' ™02 = m 22 = -^m i, m 12 = m 32 = -^ (1+4ia)a , 



m-21 = TI31 



mi2, = W33 = , TO13 and m 2 3 are O ( pr, <5 2 ) and are very small compared to the others. 



Here Diffusion constant D = 2 v f t > dimensionless Rashba coupling strength S = Afc^r and A = tpT. 
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One has to now sum four infinite series of matrices, one for each value of /?, to find out the matrix > 



using (30) iteratively. We show the explicit calculation of Y^jLo 771 i.^ below, others may be obtained following the 
same procedure. 

From Eq. (30), we find 



rn 



(2) 



= 7 
= 7 



/ 



^(? A (k', 0)m« (q, W )^(k' + q, lo) 
dk' 



(2tt)2 

it— 

7 7 (2tt)2 



mio^oaf + mii^ai^f + mu^f + m 13 g A a 3 g? 

m w g A g?a + m 11 g A g^a 1 + m 12 g A g*a 2 + m X3 g A g^a 3 



= mi raj V + mnmg Vi + ra\ 2 m\ 'a 2 + mumj V 3 
Putting mj^ in (33) and collecting the coefficients of each Pauli matrix, we obtain 

mf ] = <r [(mio) 2 + TOnTOoi + m 12 m 32 + m 13 m 23 ] + a\ [mwm w + mnm o + imu^] 

+o- 2 [mi 2 mio + mi2m 30 - «toi 3 to 2 i] + <r 3 [mi 3 mio - imiim 02 + im 12 m 31 + m 13 m 20 \ 
w (T [(»iio) 2 ] + o-i [mn(mio + to o)] + <r 2 [77112(77110 + m 30 )] + 0-3 [«mi 2 m 3 i] , 



(33) 



(34) 



where the smaller terms are neglected in the last expression. Note that, each term in coefficient of <r 3 is 0(q 2 /kp), 
but the one which is kept has much larger coefficient than the others. Similarly, 



m 



(3) 



oo(raio) 3 + o-i mn [(ttiio) 2 + (m 00 ) 2 + miomoo] 
+o- 2 [(mio) 2 + (m 30 ) 2 + m 10 m 3 o] + a 3 [im 12 m 31 (m w + m 20 + m 30 )] 

One can anticipate the higher orders now. Summing all the orders, 



(35) 



77111 



3=0 



+ <J 2 



1 - mio (1 - mio)(l - moo) 
mi2 



+ <?3 



imi 2 m 3 i 



(1 - mio)(l - m 30 ) ' (1 - mio)(l - m 2 o)(l - m 30 ) 



We similarly find 



E™o 

3=0 



Cj) 



(T - 



E 

00 

E' 

j=0 



(T - 



0) 



(TO 



+ CT1 



+ CT1 



+ CT1 



m i 



(1 


- moo)(l 


- mio) 




m 3 i 




(1 


- m 2 o)(l 


- m 30 ) 




77121 





+ <T 2 - 



+ (T2- 



™02 



(36) 

(37) 
(38) 



+ (T2 



1 - 77120 (1 - m 2 o)(l - m 30 ) 

m 22 f „ -im 21 m 32 



' (1 - 77i2o)(l - m 00 ) ' °" 3 (1 - mi )(l - m 20 )(l - m 30 ) 
It is now straightforward to calculate Xa/3 using Eqs. (31), and (36)-(39): 



(39) 



X(q,w) = 2i/iwr 



moo 
1— moo 



mpi 

(1 -moo) (1 -mio) 

m 02 



mn 

(1— moo)(l— mio) 

mio 
1— mio 

mi 2 m 3 i 



m 2 2 

(1— moo)(l— IJ120) 

m 2 im 3 2 







■■■02 A mi 2 m 3 i 

(1— m o)(l— «i2o) (l-mio)(l-m 2 o)(l-m3o) 

Q A mi2 

(1— mio)(l— m 30 ) 



(1— mi )(l-m 2 o)( 

m 2 o 
l-m 2 o 

m 2 i 



A "'21 

(l-m 2 o)(l-m 30 ) 



j m 32 

l-m 30 ) (1— mi )(l-m 3 o) 

m 3 i 

(1— m 2 o)(l— m 30 ) 

m 3 o 



1— m 3 o 



(40) 
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The diffusive propagator T> a p is connected with charge-spin density correiation function by the relation 

Xa/3 = 2viujT{D a p - 5 af3 ) 
We thus find the inverse diffusive propagator to be 



(41) 



T>~ 



1 - m 00 -ran 
-ra i 1 - m 10 

_ mn9 rn 12 m 31 







I TO12 





-im 32 
im 31 
-i m 2 i 1 - m 30 



-m 22 

"121"132 

I-TO30 
1 - ^-20 



(42) 



Instead of calculating x(q, w) explicitly as we have done above (40), if one defines 
with 



L af3 



2 



= 2 Tr 



/• dk' „ 

7 



- (1) 



e A (k',0)a /3 a K (k / + q,w) 



and 



(43) 

(44) 
(45) 

(46) 



as in Rcf. 24, one finds the inverse diffusive propagator to be 



IT 1 





-im 32 



1 - m o -ran -m 22 
-moi 1 - toio jni23 

-m 02 -imi 3 1 - m 20 im 3 i 

!»ii2 -im 2 i 1 - m 30 



(47) 



Clearly D^ 1 7^ D 12 and D^ 1 ^ 2? 21 but the other components of 2? _1 and I? are same. While V 12 and V 21 
are very small since m\ 3 and m 23 are small due to their quadratic dependence of the small parameter S, T>^ 2 and 
T>2i are not small because these are independent of S in the limit of small 5. These are indeed the important terms 
for transverse spin diffusion and hence spin-spin Hall effect. The way of obtaining inverse diffusive propagator using 
combined expressions (43)-(46) is incorrect because the noncommutative nature of Pauli matrices is not correctly 
taken care of as we have seen from our explicit calculation of \ m Eq. (40). 



ROLE OF DRESSELHAUS SPIN-ORBIT COUPLING 



Due to the bulk inversion asymmetry, the 2DES will also have Dresselhaus spin-orbit coupling. The strength of 
this coupling is however much smaller than the Rashba coupling. This appendix contains the derivation of the dc 
spin-spin Hall conductivities in presence of both Rashba and Dresselhaus spin-orbit couplings. The Hamiltonian in 
this case is 

P 2 

H = H a + V(r); H = + X (a lPy - a 2 p x ) + (3 {a lPx - a 2 p y ) (48) 

with the eigenvalues 

k-' 

e s k = — + sX RD k; X RD (6) = \J A 2 + /3 2 + AX/3 sin 6 cos 9 (49) 

777 

where (3 is the Dresselhaus SO coupling strength. The energy spectrum now depends on modified SO coupling 
strength X R d{0). Fermi surface of both the branches are anisotropic here and the corresponding Fermi momenta arc 
k s F = -smX RD + ^k 2 F + m 2 X RD . 
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The retarded (advanced) Green's functions for the electrons in this case can be written as 

rR,A ( u a 1 V" + 5(cosx k gi -sinxkQ- 2 ) , . 

9 (k ' £) = ^ e - % ± z/2r (50) 



=± 



where, cos^k = Asin ^ cosfl and sinxk = Acos f+f sin9 . 

We proceed to calculate spin-spin Hall conductivities in the same way as in section III, but the calculation is much 
more cumbersome due to the angle dependent spin-orbit coupling strength \rd{9) in the dispersion. We however find 
that one can avoid a lot of complication in self-consistent evaluation of the vertex corrected currents, by concentrating 
on the largest order in A. As long as the conductivities are non-zero in this order, the lower orders may be ignored. 
The vertex-corrected in-plane spin current components in the dc limit are thus found to be 



, 2A XB + V { JA^W-A) 

1 2 B V^ 2 - B 2 - 1 V ' 

2 1 B ^A 2 - B 2 - 1 V 1 

with A = 1 + AS 2 + 4:8%, B = 888 D , and S D = /3k F r. 

The dc spin-spin Hall conductivities, defined in Eqs. (5) and (6) can be obtained now via retarded current-current 
correlation function (7). We find 



± A ( A 2 -B 2 -A \ VA 2 -B 2 -A 

° xy " 8il B 2 ) ^J A 2 -B 2 -l ( j 

Afl\VA^W-A 



xv 8^ \B) VA 2 - B 2 - 1 

The fact that a\ y is non-zero suggests that the spin polarization of the induced spin current in the transverse 
direction is not fully perpendicular to the applied spin current. It is, in fact, fully perpendicular when one of the 
couplings is absent. Therefore in realistic situations, where /?<A, the spin polarizations of the applied and induced 
spin currents are almost perpendicular to each other as <r\ y /a^ <~ Sd/[S(1 + 4<5 2 )] <C 1 in this limit. 



[1] Datta S and Das B 1990 Appl. Phys. Lett. 56 665 

[2] Zutic I, Fabian J and Sarma S D 2001 Rev. Mod. Phys. 76 323 

[3] Awschalom D D, Loss D and Samarth N (ed) 2002 Semiconductor Spintronics and Quantum Computation (Hiedelberg: 
Springer- Verlag) 

[4] Nitta J, Akazaki T, Takayanagi H and Enoki T 1997 Phys. Rev. Lett. 78 1335 

[5] Hu C M, Nitta J, Akazaki T, Takayanagi H, Osaka J, Pfeffer P and Zawadzki W 1999 Phys. Rev. B 60 7736 

[6] Miller J B, Zumbuhl D M,Marcus C M, Lyanda-Geller Y B, Goldhaber-Gordon D, Campman K and Gossard A C 2003 

Phys. Rev. Lett. 90 076807 
[7] Murakami S, Nagaosa N and Zhang S C 2003 Science 301 1348 

Murakami S, Nagaosa N and Zhang S C 2004 Phys. Rev. B 69 235206 
[8] Sinova J, Culcer D, Niu Q, Sinitsyn N A, Jungwirth T and MacDonald A H 2004 Phys. Rev. Lett. 92 126603 
[9] Inoue J I, Bauer G E W, and Molenkamp L W 2003 Phys. Rev. B 67 033104 
Inoue J I, Bauer G E W and Molenkamp L W 2004 Phys. Rev. B 70 R041303 
Raimondi R and Schwab P 2005 Phys. Rev. B 71 033311 
Dimitrova O V 2005 Phys. Rev. B 71 245327 
[10] Mishchenko E G, Shytov A V and Halperin B I 2004 Phys. Rev. Lett. 93 226602 
[11] Chalaev O and Loss D 2005 Phys. Rev. B 71 245318 
[12] Rashba E I 2004 Phys. Rev. B 70 201309 
Rashba E I 2003 Phys. Rev. B 68 241315 
[13] Sheng L, Sheng D N and Ting C S 2005 Phys. Rev. Lett. 94 016602 

Nikolic B K, Zarbo L P and Souma S 2005 Phys. Rev. B 72 075361 Hankiewiez E M, Molenkamp L W, Jungwirth T and 
Sinova J 2004 Phys. Rev. B 70 241301 
[14] Sheng D N, Sheng L, Weng Z Y and Haldane F D M 2005 Phys. Rev. B 72 153307 
[15] Kato Y K, Myers R C, Gossard A C and Awschalom D D 2004 Science 306 1910 
[16] Wunderlich J, Kaestner B, Sinova J and Jungwirth T 2005 Phys. Rev. Lett. 94 047204 



9 



[17] Sih V, Myers R C, Kato Y K, Lau W H, Gossard A C and Awschalom D D 2005 Nature Phys. 1 31 

[18] Nomura K, Sinova J, Sinitsyn N A and MacDonald A H 2005 Phys. Rev. B 72 165316 

[19] Nomura K, Wunderlich J, Sinova J, Kaestner B, MacDonald A H and Jungwirth T 2005 Phys. Rev. B 72 245330 

[20] Sugimoto N, Onoda S, Murakami S and Nagaosa N 2006 Phys. Rev. B 73 113305 

[21] Engel H A, Halperin B I and Rashba E I 2005 Phys. Rev. Lett. 95 166605 

[22] Tse W K and Sarma S D 2006 Phys. Rev. Lett. 96 056601 

[23] Fiederling R, Keim M, Reuscher G, Ossau W, Schmidt G, Waag A and Molenkamp L W 1999 Nature 402 787 

[24] Burkov A A, Nunez A S and MacDonald A H 2004 Phys. Rev. B 70 155308 

[25] Sun Q F and Xie X C 2005 Phys. Rev. B 72 245305 
Yang H-T and Liu C arXiv:cond-mat/0604230 



